The Swift X-ray flaring afterglow of GRB 050607 
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^ ■ ABSTRACT 
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CN , The unique capability of the Swift satellite to perform a prompt and au- 

^ ! tonomous slew to a newly detected Gamma-Ray Burst (GRB) has yielded the 

^ i discovery of interesting new properties of GRB X-ray afterglows, such as the 

l/^ ! steep early lightcurve decay and the frequent presence of flares detected up to 

^ ! a few hours after the GRB trigger. We present observations of GRB 050607, 

O ! the fourth case of a GRB discovered by Swift with flares superimposed on the 

^ I overall fading X-ray afterglow. The flares of GRB 050607 were not symmetric 

"rH i as in previously reported cases, showing a very steep rise and a shallower decay, 
similar to the Fast Rise, Exponential Decay that are frequently observed in the 

O I gamma-ray prompt emission. The brighter flare had a flux increase by a factor of 

' ~ 25, peaking for 30 seconds at a count rate of approximately 30 counts s"^, and 
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it presented hints of addition short time scale activity during the decay phase. 
There is evidence of spectral evolution during the flares. In particular, at the 
onset of the flares the observed emission was harder, with a gradual softening as 
each flare decayed. The very short time scale and the spectral variability during 
the flaring activity are indicators of possible extended periods of energy emission 
by the GRB central engine. The flares were followed by a phase of shallow decay, 
during which the forward shock was being refreshed by a long-lived central en- 
gine or by shells of lower Lorentz factors, and by a steepening after approximately 
12 ks to a decay slope considered typical of X-ray afterglows. 

Subject headings: gamma rays: bursts - X-rays: individual (GRB 050607) 

1. Introduction 

The Swift Explorer (Gehrels et al. 2004), designed to discover and study Gamma- Ray 
Bursts (GRBs), was launched on 2004 November 20. When the Burst Alert Telescope (BAT, 
Barthelmy et al. 2005) triggers on a newly discovered GRB, the Swift satellite performs an 
autonomous slew to the GRB position determined on-board, allowing prompt follow-up 
observations with the two narrow-field instruments, the X-ray Telescope (XRT, Burrows et 
al. 2005a) and the Ultraviolet/Optical Telescope (UVOT, Roming et al. 2005). 

In the standard scenario the GRB emission is caused by shocks in a relativistic expanding 
fireball (Mcszaros & Rees 1997). The internal shocks of the moving ejecta are responsible 
for the prompt gamma-ray emission, while external shocks caused by the interaction of the 
fireball with the surrounding medium produce the afterglow emission observed in the X-ray 
and lower energy bands. 

Thanks to the rapid spacecraft response to the BAT triggers, the XRT has been able 
to observe and autonomously localize GRB afterglows at early times, unveiling new and 
unexpectedly complex behaviors of the X-ray lightcurves, which frequently show multiple 
breaks with segments of distinct decay slopes (Nousek et al. 2006) and flares superimposed 
on the fading afterglow. 

In this paper, wc report on the properties of the X-ray afterglow of GRB 050607, which 
was dominated, in its early phase, by emission from X-ray flares superimposed on the overall 
decay. This was the fourth case in which X-ray flares were detected in the afterglows of 
GRBs discovered by Swift (Romano et al. 2006; Godet et al. 2005; Falcone et al. 2005). 
X-ray flares have been previously observed in few cases by BeppoSAX (XRF 011030, GRB 
011121 and GRB 011211) and were identifled as the onset of the afterglow (Galli & Piro 
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2005; Piro et al. 2005). A more detailed study of the characteristics and the origin of X-ray 
flares is now possible thanks to the unique Swift capabilities, which have greatly increased the 
number of rapid observations of the early phase of the X-ray afterglows. Continued internal 
engine activity at late times was proposed to explain the very short timcscale and spectral 
variations during the flares observed in XRF 050406 (the flrst flare detected by Swift, see 
Romano et al. 2006) and the very bright flare of GRB 050502b (Falcone et al. 2005), which 
peaked 740 s after the GRB prompt emission. The flares in the X-ray afterglow of GRB 
050607 had asymmetric shapes, with a steeper rise than previously observed flares, and a 
shallower decay, similar to the FREDs (Fast Rise, Exponential Decay) that are frequently 
observed in the GRBs prompt emission. 

In this work, we describe the XRT data analysis of the GRB afterglow, focusing in 
particular on the temporal and spectral properties of the early phase of flaring activity of 
the X-ray emission. We use the notation F{t,v) oc {t — to)~°'u~^ for the temporal and 
spectral dependence of the flux. All uncertainties are given at the 90% confldence level for 
one interesting parameter (Ax^=2.7). The observations are presented in § 2 and the data 
analysis in § 3. The discussion and conclusion are in § 4. A detailed presentation of the 
procedure applied to correct for the 'pile-up' effect, caused by the high X-ray afterglow count 
rate during the flares is given in the appendix. 



2. Observations 

The Burst Alert Telescope triggered on the GRB 050607 at 09:11:22.81 UT on 2005 
June 07 (Retter et al. 2005). The spacecraft started the slew to the GRB 16.5 s after the 
trigger, while the burst was still ongoing, and settled on the burst 84 s after the trigger. 

The XRT follows an autonomous sequence of observations in response to a GRB trig- 
ger, starting with the acquisition of a short image to localize the GRB afterglow and then 
performing observations in different readout modes according to the source intensity (Hill et 
al. 2004). Due to engineering tests on the narrow fleld instruments at the time of the trigger, 
the XRT on-board centroiding algorithm could not converge on a source when the flrst im- 
age was acquired and the nominal observing sequence could not be effectively followed. The 
instrument collected a single frame in Windowed Timing (WT) mode and then immediately 
switched into Photon Counting (PC) mode, designed to provide two dimensional spatial 
information at the expense of a lower timing resolution. The GRB afterglow was monitored 
by the XRT for 30 days after the burst for a total exposure time of 363 ks. The complete 
set of the GRB observations is reported in Table 1. 
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The UVOT was in an engineering mode at the time of the GRB event and therefore 
immediate data were not available. The first UVOT observation of the GRB field started 
at 22:36:23 UT, 12 hours and 25 minutes after the BAT trigger. No optical counterpart was 
detected (Ivanushkina et al. 2005). 

The GRB 050607 field was observed with ground-based telescopes in the optical and 
the IR bands. A fading source was detected with the 4m Mayall Telescope at Kitt Peak 
Observatory in the / band with an initial magnitude of 21.46 ± 0.3, 640 s after the burst 
trigger (Rhoads 2005a). The transient was observed to have faded by 2.3±0.2 mag in the 
/ band between 640 s and 93.2 ks (Rhoads 2005c). Observations were also performed in the 
B, R, z and g bands. The afterglow was detected in the R band, but not clearly visible in 
the B band images (Rhoads 2005b), indicating a possible Lyman break in the spectrum. 

The MDM 2.4m telescope also reported a source detected in the R band (22. 5< R <23 
mag, Halpern et al. 2005a) with the observation beginning 16 min after the GRB trigger. 
The object appeared to have faded (upper hmit of i? >24 mag) in a follow-up observation 
performed on 2005 June 8 at 08:15 UT, 23 hours after the trigger (Halpern et al. 2005b). 

Observations were also performed in the IR by the REM robotic telescope (Testa et al. 
2005); prompt observations (192 s after the trigger) did not reveal any new bright object in 
the H band to a limiting magnitude of 16.0. 

3. Data Analysis 

3.1. Gamma-ray Analysis 

The BAT data (analyzed with BAT HEAsoft software v6.0.2) showed a 3-peak hghtcurve 
with a Tgo of 26.5 s, close to the mean value for long bursts of the BATSE catalog (Paciesas 
et al. 1999). The measured BAT flux had a peak of 1.15 ± 0.16 photon cm~^ s~^ and a total 
fluence (pre-slew and slew emission) of 5.8 ± 0.5 x 10~^ ergs cm~^ in the 15-150 keV band. 

The BAT spectrum of GRB 050607 is well fitted with a single power law model with 
index 1.83 ± 0.14. A fit of the spectrum with a Band function (Band et al. 1993) ^ does 
not provide a significant improvement in the xled because of the difficulties in constraining 
Epeak: duc to the narrow BAT energy range and the relative low fluence of the GRB 050607 
prompt emission. 



^f{E) = K^^{E/iQ)'^eM~E{2 + a)/Ep,,i,), if E < {a ^ f3)Ep^^^/{2 + a) and f{E) = i^3o{(a - 
/3)£;peak/[30(2 + a)]Y-^ exp(/? - a)(^/30)'', iiE>{a- /3)^peak/(2 + a). 



- 5 - 



The extrapolation of the BAT hghtcurve into the XRT energy band is performed using 
the best fit parameters obtained from the prompt emission spectral fit. If one assumes the 
GRB prompt emission to have a 'universal' spectral shape (Preece et al. 2000) described by 
a Band function with q;=1.0 and (3=2.5, the photon index of 1.83 of the simple power law 
model fit of the BAT spectrum is very steep and would yield too high a value of the flux 
in the XRT energy band. The prompt flux in the XRT band can therefore be assumed to 
lie between the fluxes extrapolated from the two spectral models: a simple power-law model 
and a Band function with fixed q;=1.0 and /3=2.5. 

The BAT spectrum was analyzed separately over different time intervals relative to the 
three peaks of the lightcurve prompt emission. In the time interval To+30 s to Tq+IOO s 
there is no source detection in the BAT observation and a 3 cr upper limit was calculated 
using the two spectral models. The BAT extrapolated fluxes are reported in Table 2. 

3.2. GRB position 

The X-ray afterglow position, determined from ground-processed data, is RA(J2000)=20'^00™42!7, 
Dec(J2000)=-F09°08'31'.'4 , with an uncertainty of 3.8". The position obtained with the xrt- 
CENTROID task was corrected for the systematic offsets observed in XRT and optical coun- 
terparts, caused by the slight misalignment between the XRT optical axis and the spacecraft 
star tracker boresight (Moretti et al. 2006). The X-ray afterglow is 44" from the reported 
BAT position (Retter et al. 2005) and 0.9" from the optical counterpart reported by Rhoads 
(2005a). The image of the afterglow detected by the XRT during the first orbit is presented 
in Figure 1, along with the BAT refined position (Retter et al. 2005), the XRT error circle 
and the optical counterpart position (Rhoads 2005a). 

3.3. Timing Analysis 

The X-ray afterglow of GRB 050607 had a complex behavior, dominated by flares in the 
early observations, followed by a phase of slow decay that lasted approximately 12 ks and a 
late steepening to values typical of GRB afterglows. The X-ray lightcurve has been extracted 
in the energy range 0.2-10.0 keV and binned to obtain 20 background-subtracted counts per 
bin. Count rate to flux conversion factors were obtained from the best spectral model flt 
for the different afterglow phases (the detailed spectral analysis is discussed in § 3.4). The 
X-ray hghtcurve of GRB 050607 and the extrapolated BAT fluxes in the 0.2-10.0 keV band 
are shown in Figure 2. 
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Due to the high count rate during the flares, the early XRT observations were severely 
affected by 'pile-up', the result of two or more photons being collected in a single CCD cell in 
a readout frame. In Appendix A we describe how the observations are altered by the pile-up 
effect and we determine the affected region and the coefficient applied to the observed fluxes 
to correct for the pile-up. 

The phase of flaring activity of the GRB 050607 X-ray afterglow is shown in Fig. 3. To 
avoid the region contaminated by pile-up, the lightcurve of the first 510 s of observations 
after the GRB trigger was extracted from an annulus that excluded events in the central 
8 pixels (18.88") and with an outer radius of 40 pixels. The first fiare (referred to as flare 'A') 
was observed from 96 s to 255 s after the GRB trigger, peaking at 4.6x10"^'^ erg cm~^s~"^ 
(8.1 counts s~^) approximately 135 s after the trigger. The second flare (flare 'B') was 
the brightest, lasting from 255 s to approximately 510 s after the trigger and including the 
majority of the X-ray emission, with a peak flux of 2.0x10"^ erg cm~^s~^ (35.6 counts s~^) ~ 
310 s after the trigger and a total fluence of 1.36x10"^ erg cm""^ (~23% of the BAT prompt 
emission) . 

The flares were asymmetric, similar to the FREDs that are often observed in the prompt 
GRB emission, with a rising side steeper than the following decay. In the case of flare A, the 
peak was reached in ~ 30 s, with the decay taking approximately 130 s, and for flare B the 
rising phase lasted for ~ 35 s, followed by a 'plateau' at the flare maximum with the X-ray 
emission at a countrate of ~ 30 counts s~^ for approximately 30 s, and a decaying phase 
~190 s long that presents marginal indications of spikes and 'bumps'. 

The flare's rise and decay can be fitted with a simple power law. With the notation 
(t — to)"i fixing to at the onset of the fiares, the fiare A slope indices are Q;A,rise=l-2±0.2 and 
ttyi,(iecay=-l-3±0.3, and for fiare B are as^rise=4.1±0.8 and Q;B_decaj/=-2.3±0.3. The choice 
of to is critical in determining the slopes indices. In fact, if to is fixed at the GRB trigger 
time ts the slopes are steeper, yielding higher values than for previously detected fiares 
(Romano et al. 2006; Falcone et al. 2005). In this latter case flare A rises as (t — ts)^'^^^'^ 
and decays as (t — tB)~^'^^°'^, while for flare B the rise and decay are fltted with (t — t^)^^^^ 
and (t - tB)-' '^° '. 

The brighter flare B decay was fltted excluding the data points of the secondary bump 
superimposed on the flare and centered 355 s after the burst trigger. From the power law flt 
decay parameters, the expected number of events in that time period is 323±18, while the 
observed number of events is 398±20. The observed 75 excess photons are signiflcant at a 
2.8 a level. Although hints of short time scale structure in the XRT lightcurve are present, 
in particular during the brightest flare, no significant periodicity was found from the timing 
analysis in the 2 — 300 s period range, using a frequency spacing of 10~^ s~^. 
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After the flares, the X-ray afterglow emission decayed to intensities for which the effect 
of pile-up was negligible. Data were extracted using a circle of 20 pixels radius (smaller than 
the one used for the flares to avoid excessive background contamination) centered on the 
GRB afterglow and the background was evaluated using a circle with r=40 pixels in an area 
of the CCD without presence of other sources. In the second, third and fourth observations 
(segments 1, 2 and 3 respectively, see Table 1) a binning of 25 events per bin was used. At 
later times multiple observations (segments) were co-added to detect the X-ray afterglow, 
and a circle of 4 pixels radius centered on the GRB position was used to obtain a higher 
signal-to-noise ratio. The measured afterglow flux was PSF-corrected to account for the 
smaller extraction region. The last two data points of the hghtcurve are 90% upper limits 
(Fig. 2). 

After the main flare, starting from 510 s after the GRB trigger, the afterglow is best 
fitted with a broken power law (shown in Fig. 2) with indices of 0.58 ±0.07 and 1.17 ±0.07, 
with a break time of 11.7±0.6 ks {Xred — 1-02). The broken power law model is a significantly 
better fit (F-test probability of 4.8x10"^) to the late X-ray afterglow than a single power 
law (index of 0.88 ± 0.03 with a xled of 1-91). The observed decay slopes and the degree of 
steepening after the break are close to the average values obtained from the sample of 27 
GRBs discovered by Swift lepoited in Nousek et al. (2006). 



3.4. Spectral Analysis 

The X-ray afterglow of GRB 050607 can be divided into four time intervals: fiare A 
extending to 255 s after the trigger, the stronger fiare B from 255 s to approximately 510 s, 
the shallow section of the afterglow emission and the late steepening after 12 ks. The spectral 
analysis (using xspec vl 1.3.2) was performed separately on each lightcurve phase, selecting 
events with grades to 12. Because of low statistics below 0.5 keV and at high energies 
the spectra were fitted in the energy range 0.5-6.0 keV; the swxpc0tol2_20010101v007.rmf 
response matrix file from the Swift XRT calibration database was used for the spectral 
fit. The ancillary file was created using the xrtmkarf command of the xrtpipeline 
software, applying the PSF correction. The spectra during the first 510 s of observation were 
extracted from an annular region as in the timing analysis to avoid pile-up contamination. 
The background was evaluated from a circle of r=40 pixels in a region without sources. 

Due to the low number of events during the first flare and the late part of the lightcurve, 
during these time intervals C-statistics (Cash 1979) were used in the fltting procedure, while 
statistics were used for the main flare. Since the X-ray count rate at low energy was too 
low to properly constrain the X-ray absorption column, the Nh was fixed to 1.41x10^^ cm~^. 
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the weighted average value in the GRB 050607 direction (Dickey & Lockman 1990). The 
flares were fitted with both a simple absorbed power law and a cutoff power law model 
to investigate a possible common origin of the X-ray flares and the prompt gamma-ray 
emission: in the energy band of the XRT data a cutoff power law is almost equivalent to a 
Band function. 

The best fit of flare A with an absorbed power law model yielded a photon index of 
2.00+ 

ai8' with the absorption fixed at the Galactic value. The brighter flare B fit with an 
absorbed power law model yielded an absorption column density 2a higher than the Galactic 
value. To allow comparisons with earlier and later spectra, the absorption was fixed to the 
Galactic value, obtaining a photon index of 2.27+q;J2) while selecting only the decaying phase 
of the flare yielded a photon index of 2.311q ^g. 

The fit of flare B with a cutoff power law model yielded a best fit absorption column 
density closer to the Galactic value than in the single power law case. Fixing the absorption 
to the Galactic value, the cutoff power law model is a better fit than an absorbed power law 
{xled improved from 1.54 to 1.14). The parameters of the cutoff power law model fit are not 
very well constrained, especially in the case of the cutoff energy, due to the limited number 
of photons in the dataset at high energies. The cutoff power law model fit yielded a photon 
index of 1.411q'53, harder than the simple power law, and a cutoff energy of 2Atll keV. 
Our attempt to also fit the dimmer, first fiare with a cutoff power law did not properly 
constrain the fit parameters due to the low number of photons during that time interval. 
However, we investigated the spectral evolution of the afterglow during the fiarcs, dividing 
the lightcurve into a low energy band and a high energy band and evaluating the band 
ratio. The energy ranges (0.2-1.5 and 1.5-10.0 keV for the low and the high energy bands 
respectively) were selected to obtain a similar number of counts in the two bands. Spectral 
evolution is evident in the band ratio plot (see Figure 4). In particular, the hard band 
contributed to the overall emission mostly during the rising part of the two flares, while the 
soft emission lingered longer and started decaying at later times. The observed behavior 
indicates that the emission during the fiares is harder than the underlying afterglow, and 
that the flares evolve spectrally, softening as they decay. 

The late part of the X-ray afterglow was analyzed using events from a circular region 
of 20 pixels radius. We applied C-statistics to flt the data with an absorbed power law 
with absorption fixed to the Galactic value, yielding a photon index of l.TS+o'lg for the 
shallow section of the lightcurve and a photon index of 1.97lg 34 after the lightcurve break. 
The spectra of the different afterglow phases are shown in Figure 5 and the fit results are 
presented in Table 3. 
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4. Discussion and Conclusions 

The X-ray lightcurvc of GRB 050607 has a complex behavior, presenting most of the 
unexpected features that have been observed in the afterglows of GRBs discovered by Swift, 
such as the presence of multiple breaks in the lightcurve and the flaring activity superimposed 
on the overall lightcurve decay. 

XRT observations indicate that flares are common in the early phase of X-ray after- 
glows (Burrows et al. 2005b; Zhang et al. 2005; Nousek et al. 2006): as of 2005 October 20, 
flares were detected by the XRT in 23 of the 50 afterglows for which observation with the 
narrow field instruments started within 10 minutes from the BAT trigger. Several scenarios 
have been proposed as the possible source of fiux variabilities of the X-ray afterglows: den- 
sity fiuctuations in the enviroment (Lazzati et al. 2002); angular structures in the outfiow 
(Nakar et al. 2003); refreshed decelerated shocks by slower shells (Rees & Meszaros 1998); 
synchrotron emission from a reverse shock (Kobayashi et al. 2005). The time variability 
for these models, all related to external shock emission, is longer than the observed time 
since the GRB trigger (loka et al. 2005). Therefore, they are unable to account for the very 
short timescales of the two flares of GRB 050607, that showed in particular a very steep 
rise, sharper than for previously reported cases (Romano et al. 2006; Falcone et al. 2005), 
with values of Strise/tpeak ^ 0.2, and slower dacays with Stdecay/tpeak ^ 1- The flares of GRB 
050607, unhke the ones previously detected by the XRT, did not have a symmetric shape 
but were rather similar to the FREDs that are common in the prompt GRB emission, and 
had indication of additional short time variability superimposed on the flares emission. 

The choice of the spectral model to flt the X-ray flares observed is important for a 
correct interpetation of the data: in the case of GRB 050607 the cutoff power law model 
(coincident to a Band function in the XRT energy band) provided a significantly better fit of 
the brighter flare spectrum than a simple absorbed power law. In addition, the hardness ratio 
showed that the emission was harder at the flare onset, while the decay was characterized by 
a softer component. In the late flares detected by BeppoSAX, the flare spectral properties 
were similar to those of the late afterglow, and they were therefore interpreted as the onset 
of the afterglow itself. The spectral evolution observed during the flares of GRB 050607 
suggests instead a mechanism for the origin of the flares distinct from the one responsible 
for the emission of the underlying decaying afterglow. 

The described characteristics observed in the flares of GRB 050607 strongly indicate 
internal shock emission from late central engine activity at a few hundred seconds after the 
prompt gamma-ray burst as the origin of the flares. Compared with the three peaks in 
the prompt emission detected by the BAT, the X-ray flares had lower energy and longer 
time scale. If the late internal shocks happened at larger radii than those of the prompt 
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gamma-ray emission, smaller magnetic field strength (internal energy density) and a larger 
curvature time scale (angular spreading time scale) for the shocks would make the resultant 
flares softer and longer, consistent with our observations. At larger radii the fireball shells 
should spread and have a larger width, and the observed plateau of the bright flare could be 
interpreted as produced by a thick shell collision or by a series of collisions. 

For a large fraction of Swift GRBs, the X-ray afterglow is characterized by a steep to 
shallow to 'normal' behavior (Chincarini et al. 2005; Tagliaferri et al. 2005; Zhang et al. 
2005; Nousek et al. 2006; Panaitescu et al. 2005). In the case of GRB 050607 the flares are 
presumably superimposed on the steep decay and the transition to the shallow phase is not 
well constrained. A shallow to 'normal' transition was observed ~ 12ks after the burst. The 
steepening (q;=0.58±0.07 to 1.17±0.07) is larger than the value of Sa—l/A expected when 
the cooling frequency passes through the observed band in the standard constant energy 
blast wave model and therefore implies a change of the hydrodynamics of the forward shock, 
suggesting that the forward shock continues to be refreshed for some time (Zhang et al. 
2005; Nousek et al. 2006). In one possible mechanism, the bulk Lorentz factor of the ejecta 
decreases at the end of the prompt phase, resulting in a monotonic increase of the Lorentz 
factor with radius. This outflow gradually catches up with the forward shock (blast wave), 
and gives a smooth energy injection (Sari & Meszaros 2000; Panaitescu et al. 2005). In 
another scenario, the energy injection originates from a long-lived central engine, active up 
to a few hours after the prompt emission (Dai & Lu 1998). The decay slope after the break 
{a— 1.17) is consistent with the standard interstellar medium (ISM) forward shock model 
(Sari et al. 1998). If we apply the standard blast wave model to explain the light curve after 
12ks, an electron index p ~ 2.2 is obtained from the decay and spectral indices. 

If one assumes that the shallow decay observed in the X-ray energy band originates 
from refreshed shocks, emission from the associated reverse shocks is expected to peak in 
the optical band (Sari & Meszaros 2000), with a decay shghtly slower than the forward 
shock emission. The / band decay during the first 6400 s could be modeled with a single 
power law of index 0.5 (Rhoads 2005b). The measured / magnitudes and the X-ray fluxes 
are consistent within the uncertainties with the model predictions during the shallow phase. 
With the late measurement of the transient, still detectable in the / band 93 ks after the 
burst trigger (Rhoads 2005c), the / band afterglow decay could be modeled with a single 
power law decay with slope 0.4, or 0.3 if the first detection after 640 s was excluded. The fit 
of the X-ray afterglow decay after the flares with a single power law yields a slope of 0.88. 
The difference between the reported / band decay and the X-ray slope suggests that the late 
/ band detection had a significant contribution coming from the GRB host galaxy. A new 
deep observation of the field with ground based telescopes would verify this interpretation. 
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A. Pile-up correction 

The pile-up effect causes two or more photons to be considered by the CCD as a single 
event with energy equal to the sum of the energies of the individual photons. Pile-up, 
therefore, has a major impact on the observations, by lowering the apparent source count 
rate and making the apparent spectral index harder (for a detailed discussion see for example 
Weisskopf et al. 2000). 

Pile-up also modifies the grade distribution of the photons collected on the CCD, via 
an effect known as 'grade migration' (Davis 2001). When an X-ray photon hits the CCD 
surface, the charge produced is spread over one or more pixels; single pixel events are recorded 
as 'grade 0' events, split pixel events are designated with higher grades (1 to 12 for X-ray 
events, see Burrows et al. 2005a for the XRT grade definitions) . The effect of pile-up is to 
decrease the number of single pixel events and to increase the number of apparent multiple 
pixel events, depending on the overall incident count rate. 

To investigate the extent of the CCD region affected by pile-up we derived the grade 
distribution of the events collected from the GRB afterglow during the main fiare, when the 
XRT count rate exceeded 1 count s~^. Annular regions centered on the GRB emission with 

increasing inner radii and fixed outer radii {rout = 40 pixels) were extracted to verify the 
spatial extent to which the observations are affected by pile-up. The ratio of single pixel 
events to the total number of events has been calculated for each annular region and the 
results arc reported in Table 4. For comparison, ground calibration tests (Morctti et al. 2004) 
with the XRT showed that in non piled- up regimes the fraction of single pixel events is 0.78 
at 1.49 keV (the mean photon energy during the bright fiare, once the pile-up correction was 
performed, was measured to be 1.7 keV). 

For annular regions with inner radii less than 5 pixels the 'grade migration' effect caused 
by pile-up is significant, while for inner radii larger than 8 pixels the fraction of single pixel 
events remained constant and became comparable with ground calibration measurements 
for sources with low count rates. To confirm the result, the radial profile of the bright fiare 
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was extracted and compared with the analytical PSF described in Moretti et al. 2004. The 

two profiles showed little or no deviations at radii larger than 8 pixels. To avoid the region 
affected by pile-up while minimizing the number of rejected events in the central region, the 
analysis of the GRB afterglow during the bright flare was performed using an annular region 
that excluded events in the central 8 pixels (18.88") and with an outer radius of 40 pixels. 

The number of events rejected from the central 8 pixels is in principle dependent on 
the incident spectrum and the details of the energy-dependent PSF. For a point source at 2' 
from the center of the CCD (as is the case of the GRB 050607 X- ray afterglow during the 
first orbit), 80% of the PSF is included in an 8 pixel circle at an energy of 0.5 kcV, 77% at 
1.7 keV (the mean energy of the fiare spectrum corrected for pile-up) and 74% at 3.0 kev. 
We therefore used XSPEC vl 1.3.2 to calculate the correction factor for the pile-up. The 
spectra from the annular region and the background region were extracted; the xrtmkarf 
command (version 0.4.14) was used to create two ancillary response files, with and without 
applying the PSF correction. Fitting the spectrum with an absorbed power-law model we 
obtained the best-fit parameters for the two ancillary response files. The ratio of the two 
normalization parameters yields the correction factor to account for the effect of pile-up 
. The coefficient {cpHeup = 0.20 ± 0.03) obtained from the spectral analysis was used to 
multiply the count rate observed in the annular region and is consistent with the fraction of 
PSF contained in the 8 pixel central region. 
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Fig. 1.— The GRB 050607 afterglow XRT image extracted from the first orbit of Photon 
Counting observation. The BAT refined ground position, the XRT 3.8" error circle and 
the optical afterglow position (white cross inside the XRT error circle) determined with the 
Mayall Telescope at Kitt Peak (Rhoads 2005a) are shown. 




Fig. 2.— The combined BAT and XRT GRB 050607 lightcurve. The BAT fluxes were 
extrapolated to the 0.2-10.0 keV energy band assuming a simple power law (squares) and a 
Band (asterisks) spectral model. In the last BAT time bin (30-100 s) there was no detection 
in the BAT data and upper limits were derived. The X-ray afterglow from the entire set of 
observations is in the energy band 0.2-10 keV, pile-up corrected. The time bars correspond 
to the time bin sizes. The best-fit broken power law model is shown. 
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Fig. 3. — The XRT lightcurve of the early GRB 050607 afterglow, in the energy band 
0.2-10 kev, corrected for pileup, with 20 photons per bin. The early X-ray emission was 
characterized by the presence of flares. Flare A had a peak at ~ 8 counts the brighter 
flare B lasted from 255 to 510 s after the first trigger with emission around 30 counts for 
approximately 30 s. Both flares are asymmetric: the decay is shallower than the very steep 
rise. This may be due to additional short timescale activity superimposed on the main flare, 
like the small 'bump' observed at approximately 355 s after the burst trigger. 
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Fig. 4. — The band ratio of the early phase of the X-ray afterglow. The hard (H=1.5- 
10.0 keV) and the soft (8=0.2-1.5 keV) band energy ranges were chosen to obtain a similar 
number of counts in the two bands. The band ratio (H/S) shows evidence of spectral evolu- 
tion during the two flares, with a harder emission at the flare's onset and a soft component 
that lasted longer and dominated the decay of the flare. 
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Fig. 5. — The XRT spectra of the different afterglow phases. The data were fitted in the 
energy range 0.5-6.0 keV with a binning of 20 photon per bin for the bright second flare 
and a binning of 1 photon per bin for the other phases, for which C-statistics was used. The 
data are rebinned here to 10 photons per bin for presentation purposes, and the first flare 
spectrum is normalized by a constant factor of 0.2 for clarity. The absorbed cutoff power 
law model was used to fit the spectrum of the bright fiarc while a simple absorbed power 
law model was used for the other afterglow phases. The cutoff power law model yielded a 
harder spectrum (r=1.41) for the bright flare than the late afterglow. In the lower panel 
the differences between the data and the fitted models are displayed. 
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Table 1. Observations of GRB 050607 



Sequence Number 


Number of orbits 


Start Time 


Stop Time 


Exposure Time 






UT 


UT 


N 





3 


2005-158-09:12:56 


2005-158-13:04:40 


7955 


1 


13 


2005-158-17:13:52 


2005-159-13:11:19 


44423 


2 


11 


2005-160-00:05:00 


2005-160-16:30:51 


30730 


3 


13 


2005-161-00:05:00 


2005-161-23:03:17 


26397 


A 


c 
o 


zuuo- iuz-uu.u^. oy 


ZUUO-lUZ-Zl.UO.O 1 




5 


3 


2005-165-19:38:15 


2005-165-23:29:47 


7721 


6 


14 


2005-167-00:36:29 


2005-167-23:42:54 


31998 


7 


12 


2005-168-02:17:51 


2005-168-23:44:58 


24316 


8 


29 


2005-169-23:17:40 


2005-171-23:59:57 


39943 


9 


14 


2005-172-01:22:51 


2005-172-23:59:59 


12775 


10 


39 


2005-176-01:34:52 


2005-178-21:41:30 


62773 


11 


12 


2005-181-01:01:56 


2005-181-23:37:51 


12189 


12 


3 


2005-182-01:07:13 


2005-182-12:29:26 


4160 


13 


3 


2005-182-17:15:46 


2005-182-20:31:42 


1682 


14 


11 


2005-182-23:41:26 


2005-183-23:51:23 


7676 


15 


9 


2005-184-01:23:59 


2005-184-22:21:44 


6979 


16 


21 


2005-184-23:45:48 


2005-186-12:55:57 


21884 


17 


6 


2005-186-14:25:03 


2005-186-22:35:26 


5005 


18 


11 


2005-187-00:04:59 


2005-187-22:42:18 


10566 
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Table 2. GRB 050607. BAT extrapolated Fluxes 



Time range 

N 


Power-law Flux 
[10~^° erg cm~^s~^] 




Band Function Flux 
[10-^° erg cm-2s-i] 




-3.0-10.0 
10.0-16.5 
16.5-30.0 

30.0 100.0 


260.581^^?^ 
470.141^^^:^^ 

<230 


0.89 
1.21 
0.73 


41.951^^7^ 
51.89t|?^ 

52.80_22!30 

< 0.51 


0.93 
1.27 
0.72 



Note. — GRB 050607: extrapolation into the 0.2-10.0 keV band 
of the BAT data from the prompt emission; two spectral models were 
assumed: a power-law and a Band Function. The Xred '^^ ^^c spectral 
model best-fit is reported. For the fourth time interval, where no source 
is detected in the BAT data, an upper limit was derived for each of the 
two spectral models. 



Table 3. GRB 050607. X-ray afterglow: spectral and timing parameters 



Phase Time from trigger T xfed ^ Xred P " (decay slope) 

Power Law Cutoff Power Law Power Law 



First flare 


96s<T<255s 


2 nn+°-i^ 

^•^^-0.18 


142.2'' 






Second flare 


255s<T<510s 


9 97+0.13 
^•^'-0.12 


1.54 


1 41+0.48 -.14 

-'-•^-'^-0.53 ^-^^ 




Second flare decay 


315s<T<510s 


9 O1+0.21 

^•^■■--o.ig 


0.89 


1 01+0.21 
^•''^-o.ig 


7.20 (2.26 if To = Tfla.e) 


Flat 


510s<T<12ks 


1 70+O.I8 
■'-'*-0.13 


150.5'' 


0.78l°:i« 




Break 


T>12ks 


1 Q7+0.36 
-•-•^'-0.34 


127.8^ 


• • • 0.971°;^ 


1 17+0.07 



^Value obtained with C-statistics fit 

Note. — Spectral and timing parameters of the X-ray afterglow evolving phases. The spectra were extracted in the energy 
range 0.5-6.0 keV with an absorbed power law fit. The Nh was fixed to the Galactic value of 0.14xl0^^cm~^ due to insuflacient 
statistics at low energies to properly constrain the X-ray absorbing column. We investigated the possibility of internal shocks (the 
mechanism responsible for the prompt emission) as the origin of the bright second flare, modelling the spectrum with an absorbed 
cutoff power law and obtaining a better fit to the data than using a simple power law. As can be seen from the photon indices, 
the choice of the model is critical for determining the hardness of the spectrum. Not enough counts were collected during the first 
flare to constrain the absorbed cutoff power law parameters in that time interval, and a simple absorbed power law was used. 
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Table 4. GRB 050607. Grade distribution 



Inner Radins 


E\Tnts 


Single ]:)ixel events 


Fraction of total e\Tnts 


pixel 


1774 


1131 


63.7% 


5 pixels 


1092 


805 


73.7% 


6 pixels 


958 


726 


75.8% 


7 pixels 


806 


631 


78.3% 


8 pixels 


718 


569 


79.2% 


9 pixels 


634 


511 


80.6% 


10 pixels 


567 


456 


80.4% 


15 pixels 


362 


292 


80.6% 



